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ABSTRACT 

We present high-resolution spectroscopy of the oxygen if-shell interstellar absorption edge in seven 
X-ray binaries using the High Energy Transmission Grating Spectrometer (HETGS) onboard the 
Chandra X-ray Observatory. Using the brightest sources as templates, we found a best-fit model of 
two absorption edges and five Gaussian absorption lines. All of these features can be explained by the 
recent predictions of _fT-shell absorption from neutral and ionized atomic oxygen. We identify the Ka 
and Kj3 absorption lines from neutral oxygen, as well as the 5 = 3/2 absorption edge. The expected 
5 = 1/2 edge is not detected in these data due to overlap with instrumental features. We also identify 
the Ka absorption lines from singly and doubly ionized oxygen. The O I Ka absorption line is used 
as a benchmark with which to adjust the absolute wavelength scale for theoretical predictions of the 
absorption cross-sections. We find that shifts of 30-50 mA are required, consistent with differences 
previously noticed from comparisons of the theory with laboratory measurements. Significant oxygen 
features from dust or molecular components, as suggested in previous studies, are not required by 
our HETGS spectra. With these spectra, we can begin to measure the large-scale properties of 
the interstellar medium (ISM). We place a limit on the velocity dispersion of the neutral lines of 
<200 km s , consistent with measurements at other wavelengths. We also make the first measurement 
of the oxygen ionization fractions in the ISM. We constrain the interstellar ratio of O II/O I to «0.1 
and the ratio of O III/O I to <0.1. This work demonstrates the utility of X-ray spectroscopy for 
studies of the ISM. Future work will provide measurements of the relative abundances and ionization 
fractions for elements from carbon to iron. 

Subject headings: ISM: general — X-rays: ISM — X-rays: binaries 



1. INTRODUCTION 

In the space between the stars resides the interstellar 
medium (ISM), composed of gas and dust and containing 
both compact clouds and a diffuse component. The dif- 
fuse ISM has four major phases: cold neutral, warm neu- 
tral, warm ionized, and hot ionized (Heiles & Ku lkarnJ 
Il987|) . Despite years of study, questions still remain as 
to the relationship between these components. High- 
resolution X-ray observations from the Chandra X-ray 
Observatory and XMM-Newton offer a new window with 
which to study the ISM. 

The ISM affects X-ray spectra in two ways: photoelec- 
tric absorption, particularly at low energies, and scat- 
tering by dust grains, producing X-ray halos. At X-ray 
energies, absorption features are from the excitation and 
ionization of inner-shell (AT-shell) electrons, although for 
high-Z elements like iron, L-shell absorption edges are 
also detectable. The wavelength range available to Chan- 
dra and XMM includes absorption features from carbon 
through iron. Since the abundances of these elements arc 
substantially lower than for hydrogen and helium, X-ray 
studies can probe the ISM over larger distances than is 
possible in the optical and ultraviolet. Measurements of 
the ISM absorption features in the spectra of X-ray bina- 
ries will allow us to study the abundances and ionization 
fractions for a large number of elements. This study is 
analogous to optical and ultraviolet observations of stars 
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and quasars that measure ISM and intergalactic medium 
absorption features. X-ray binaries are the brightest ob- 
jects in the X-ray sky and therefore require only mod- 
erate observation lengths to produce high signal-to-noise 
spectra. 

An understanding of the ISM absorption features 
is also important for determining the intrinsic spec- 
tral properties of X-ray sources. Early photoelectric 
absorption models assumed that the ISM was com- 
posed of a neutral, monoatomic, h omogeneous gas 
with s olar abundance s (se e e.g., iBrown fc Gouldl 
I1970D . iRide fc Walker! l|1977ft were the first to com- 
pute a more complex 2-phase photoelectric absorption 
model which more closely resembles the actual ISM. 
Improvements in X-ray absorption cross-sections, 
and solar and ISM abund ance data, led to up- 
dated ab sorption models |M orrison fc McCammonl 

^83{ |Balucinska-Ch urch fc McTjtmimonl Il992t 

IWilms. Allen, fc McCravll200fJ() . 

While the recent models provide good fits to low- 
and medium-resolution X-ray spectra, they are not 
appropriate for the grating spectra available with 
Chandra and XMM. Most importantly, the absorp- 
tion cross-sections utilized in t h e standard models 
llHenke. Gullikson. fc Da"jJisl Il993t iVerner et all IT993T: 
Verncr & Yak ovlevl I1995D do not accurately describe 
the resonance structure of the edges. Instead, the 
cross-sections are modeled by a step function at a 
fixed edge position. The positions were d etermined ei- 
ther from low-resolution solid-state data ijHenke et alJ 



2 



JUETT ET AL. 



19931) or from theoretical calculations ijVerner et al.ll99"3l 



Verner fc Yakovlevl 1 19951k neither of which can be reli- 
ably taken as the true atomic edge position. Addition- 
ally, the models do not include any contribution from ion- 
ized forms of the elements. The limitation s of the stan- 
dard absorption models are well known fsee lWilms et al.l 
2000) , yet no high-resolution model is available for stan- 
dard usage. 

The first attempt to measure ISM absorption edges in 
the X-ray s used the Einstein Focal Plane Crystal Spec- 
trometer (Schattcnburg fe Canizaresll986|) and found ev- 
idence for the O I ls-2p (Ka) absorption resonance and 
a possible O II edge. Recent Chandra and XMM results 
have revealed more complex structure around the oxygen 
i^-shell absorption edge. Three studies have used Chan- 
dra to study ISM absorption edges in X-ray binary spec- 
tra llPaerels et al.ll200lt iSchulz et all l2002t ITakei et all 
2002). Each found a prominent O I ls-2p resonance ab- 
sorption line at 23.51 A and evidence for another line 
feature at «23.36 A, attributed to compound forms of 
oxygen, but also co nsistent with the O II ls-2p transition. 
ITakei et"aTI (|2002f) found that the oxygen edge in Cyg X- 
2 could be best fit by a model consisting of three absorp- 
tion edges, representing the two edges in atomic oxygen 
absorption and a contribution from oxygen in molecular 
forms. Using the X MM Reflection Grating Spectrometer, 
Ide Vries et all l|2003l) found that the ISM oxygen K-sheW 
edge of both X-ray binaries and extragalactic sources was 
well described by the theoretical R-matrix calculation 
of atomic oxygen l|McLaughlin fc Kirbvl lT998?h although 
requiring a shift of the edge complex relative to the ls- 
2p transition. All of the high-resolution work points to 
a complex system of absorption features at the oxygen 
edge, not the simple step function used in absorption 
models. 

The purpose of this work is to compile a large sam- 
ple of spectra to study the oxygen K-sheW edge complex. 
Other ISM edges will be studied in a future work. One of 
our goals is to provide the community with a reasonable 
model for the oxygen if-shell edge complex for use in 
future analyses. The high-resolution astrophysical data 
can also be used to test our understanding of the atomic 
structure that leads to features in photoelectric absorp- 
tion by comparing with theoretical calculations. Finally, 
we will study the nature of the ISM by measuring the 
abundances in various ionizations states and studying the 
uniformity of the ISM absorption. Combined with simi- 
lar studies of other edges, high-resolution X-ray spectra 
will in the future probe relative abundances and deple- 
tion in the ISM. 

2. EXPECTED STRUCTURE OF THE OXYGEN 
ABSORPTION EDGE 

The interstellar absorption features found in the X-ray 
band are due to ionization of inner-shell (n = 1) electrons 
in neutral atoms. This absorption results in a series of 
absorption edge complexes, one for each element. The 
cross-section for each edge has a similar large-scale en- 
ergy dependence. Photons with energies below the ion- 
ization potential cannot ionize the atom, and the cross- 
section is negligible. At the ionization energy, the cross- 
section is large and then decreases roughly as E~ 2 - 7 at 
higher energies. The reduction at higher energies is due 
to the increasing probability that a photon has enough 
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Fi g. 1. — Calculated theoretical cr oss-section of neutral oxygen, 
from Gorczvca & McLaughlin (2000). Labeled are the positions of 
the ls-2p and ls-3p transitions and the edges. The dotted lines 
show the positions of the 5 = 3/2 and S = 1/2 resonance series 
limits at 22.781 and 22.576 A, respectively. 



energy that it just passes by the atom without ionizing it. 
This description is employed in the standard absorption 
models and is reasonable for CCD-resolution spectra. 

But with the higher-resolution spectra available from 
grating instruments, we can resolve finer structure. The 
most important effect is due to excitation of a X-shell 
electron into a higher orbital. Photons with energy 
slightly less (or wavelengths greater) than the ioniza- 
tion energy can excite the electron into a higher n shell. 
This produces narrow absorption lines at the excitation 
energies yielding the \s-2p (Ka), ls-3p (K/3), etc., fea- 
tures. The ls-np series of lines is termed a resonance 
series, ending at the series limit or ionization energy 
(n — > oo). More structure in the edge results from con- 
sidering the spin of the final state of the atom. For oxy- 
gen, the ionized + product with a Is hole can have 
total spin angular momentum S — 1/2 or S — 3/2, 
producing separate series limits at different energies due 
to the spin-dependent exchange force. Each series limit 
has an associated resonance series from neutral oxygen 
excitations which have the same configuration of the 
inner shells (n = 1,2 with a Is hole), neglecting the 
contribution of the valence (n > 3) electron. In ad- 
dition, differences in the total orbital angular momen- 
tum L of the final state also give rise to different reso- 
nance series. The dominant final states have L = 1 with 
5=1/2 and S = 3/2, the configurations ls2s 2 2p 4 ( 2 P) 
and ls2s 2 2p 4 ( 4 P), respectively. There are also addi- 
tional series associated with the L — and L = 2 fi- 
nal states, ls2s 2 2p 4 ( 2 5*) and ls 2s 2 2jp 4 ( 2 £>), respectively 
fsee iMcLaughlin fc Kirbvlll998|) . but these will produce 
weaker features since they require both a promotion of a 
ls-np electron and a transition in the 2p 4 subshell. These 
weaker features are undetectable with current X-ray ob- 
servatories. 

To calculate the theoretical photoelectric absorption 
cross-sections, one must solve the Schrodinger equation. 
For atoms with only one electron (hydrogen-like) this cal- 
culation is straightforward and exactly solvable, but as 
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Table 1. Theoretical Predictions for Oxygen K Edge Features 





Predicted 


Shifted 


Feature 


A (A) 


A (A) 


Gorczyca & 


McLaughlin 


2000 


ls-2p 


23.532 


23.508 a 


ls-3p 


22.907 


22.884 


ls2s 2 2p 4 ( 4 P) 


22.781 


22.758 


ls2s 2 2p 4 ( 2 P) 


22.576 


22.553 


Pradhan et al. 2003 


O 1 \s-2p 


23.45 


23.508 a 


O II ls-2p 


23.27 


23.33 


O III ls-2p 


23.08 


23.14 




23.02 


23.08 




22.93 


22.99 


O IV l«-2p 


22.73 


22.79 




22.67 


22.73 




22.46 


22.52 



a Referenced to weighted mean of observational results, see § 5. 



the number of electrons in the atom increases, so does the 
complexity of the calculation. Multi-electron atoms re- 
quire numerical methods to solve for the cross-sections. 
Such methods are only able to approximate the solu- 
tion of the Schrodinger equation. T here are a num ber 
of numerical techniques employed fsee lBau tista 20Qjj for 
a review). The preferred method for calculating cross- 
sections that include resonance structure is the R-matrix 
method. While the structure and strength of the cross- 
sections calculated using the R-matrix method compare 
well with laboratory measurements, the absolute posi- 
tions are different. This results from the truncation of 
the wavefunction description of the resulting ion, which 
is necessary to keep the computation manageable. In 
practice, the energies of the theoretical cross-sections are 
shifted when compared to laboratory data. 

Motivated by the high-resolution spectroscopy pos- 
sible with Chandra and XMM, several groups have 
published theoretical calculations of the cross-sections 
of inner-shell processe s in a to mic and ionized oxygen 
iMcLaughlin fc Kirbvl 119981 iGorczvca fc McLaughlin! 
120001: iPradhan et alJl2003[) . As mentioned above, the 
cross-section of atomic oxygen in the vicinity of the 
K edge is dominated by two resonance series con- 
verging to the limits cor responding to the 2 P an d 
4 P states of 0+ f see e.g.. IMcLaughlin fc Kirbvl f l998). 
Laboratory data tMeiizel et alJl!996t iStolte et al.lll997t 
IMcLaughlin fc Kirbvll 9981 compare well qualitatively to 
these calculations, but discrepancies in the positions of 
the re sonan ce features were noted by b oth IStolte et alJ 
l|1997ft and IMcLaughlin fc Kirbvl lll9 98lK More recent 
work (Gorczvca & McLaughlin 2000) has improved on 
the t heoretical calculations in order t o better match the 
data. IGorczvca fc McLaughlin! l)2000j) adjusted their the- 
or etical calculat i on to match the threshold energies found 
bv IStolte et alJ l)1997j) . producing theoretical data that 
matches laboratory measurements well in both cross- 
section and position of the resonance features. It should 
be noted that while laboratory measurements have very 
good relative accuracy, the absolute positional accu- 
racy is less well known. We therefore find it reason- 
able to suggest that the energy scale of the laboratory 
measurements, and consequently the theoretical calcula- 
tions, may require a shift to match the wavelength cal- 



ibration of the Chandr a data. Figure 1 shows the cal- 
culated cross-section of IGorczvca fc McLaughlinl ( 2000) 
along with identified features. Table 1 lists the predicted 
wavelengths of the important resonance features, i.e. the 
\s-2p and ls-Z p resonance lines, and the two se ries lim- 
its as found bv IGorczvca fc McLaughlinl l)2000l) . as well 
as the expected values of those features when the ls-2p 
transition is shifted to the mean value found in our anal- 
ysis (see §5). 

Along with neutral oxygen, we would expect i onized 
oxygen to be present in the ISM. IPradhan et al.l 1)20031) 
calculated the positions and cross-sections of the ls-2p 
resonance features from all ionized forms of oxygen. Ta- 
ble 1 lists the predicted wavelengths of the O 1-0 IV 
lines, as well as their positions when the O I ls-2p po- 
sition is shifted t o the mean value fou nd in our anal- 
ysis (see §5). The IPradhan et all II2003D result shows a 
signifi cant discrepancy with the lGorczvca fc McLaughlinl 
( 2000) result for the predicted wavelength of the O I ls- 
2p transition. While the absolute positions may differ 
due to differences in the calculation, we assume that the 
relative positions of these features are reliable. We can 
test this assumption by comparing the theoretical posi- 
tion of t he O II line to the lab oratory measurement of 
23.35 A ifkawatsura et alJl2002|) . The position given by 
IPradhan etafl l|2003ft is 23.27 A, but when our shift is ap- 
plied, the theoretical position becomes 23.33 A, in much 
better agreement with the laboratory measurement. We 
have recently become aware of new calculations of the full 
iv'-shell absorption edge cross-sections for ionized oxygen 
(Gorczyca 2004, in prep.), but do not include these re- 
sults in our analysis. 

3. OBSERVATIONS AND THE HETGS OXYGEN EDGE 

Our study is based on archival observations of seven 
bright X-ray binaries obtained with Chandra using 
the High Energy Transmission Grating Spectrometer 
(HETGS) in combination with the Advanced CCD Imag- 
ing Spectrometer (ACIS; l^iitoifetd] E5^). The 
HETGS spectra are imaged by ACIS, an array of six 
CCD detectors. The HETGS/ ACIS combination pro- 
vides both an undispersed (zeroth order) image and dis- 
persed spectra from the gratings. The various orders 
overlap and are sorted using the intrinsic energy resolu- 
tion of the ACIS CCDs. The HETGS carries two trans- 
mission gratings: the Medium Energy Gratings (MEGs) 
with a range of 2.5-31 A (0.4-5.0 keV) and the High 
Energy Gratings (HEGs) with a range of 1.2-15 A (0.8- 
10.0 keV). The first-order MEG (HEG) spectrum has a 
spectral resolution of AA = 0.023 A (0.012 A). 

Of the p revious ISM studies of the oxygen if-shell edge, 
only one (Schulz ct al. 2002J) was performed with the 
HETGS. The oxygen edge is accessible with the HETGS 
exclusively via the MEG. The other studies were per- 
formed with either the Low Energy Transmission Grating 
Spectrometer (LETGS) on Chandra, or the Reflection 
Grating Spectrometer (RGS) on XMM. The MEG of- 
fers the highest available spectral resolution A/AA«1000, 
compared with 460 for the LETGS, and 300 for the RGS 
at 23 A, but the lowest effective area (5, 8, and 60 cm 2 
for the MEG, LETGS, and RGS respectively at 23 A). 

The MEG images both the positive and negative grat- 
ing orders onto ACIS. For the oxygen edge region of the 
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Fig. 2. — Effective area curves of the MEG —1 order spectra on 
ACIS-S1 for two observations. The solid line is from the Cyg X-2 
observation taken 1999 September 23, while the dashed line is from 
the Cyg X-l observation taken 2002 July 30 (ObsID 3724). The 
difference in the effective area is due to the contaminant which is 
building up on the ACIS detector. 



spectrum, the —1 order is readout by the back-side illu- 
minated SI CCD, while the +1 order is readout by the 
front-side illuminated S4 CCD. In the brightest sources, 
we found a count excess in the +1 spectrum at the oxygen 
edge, possibly from a previously undetected hot pixel, 
and have therefore chosen to restrict our analysis to the 
— 1 MEG spectrum. Since the effective area of SI is «6 
times larger than S4, this does not significantly affect our 
countrate. The effective area curve of the MEG —1 order 
is shown in Figure 2. The structure in the oxygen if -shell 
region is from internal absorption in the CCDs and the 
optical blocking filter (OBF). The feature at 23.33 A is 
likely from the ls-2p resonance of double-bonded oxy- 
gen found in the OBF, which is composed of polyimidc 
(a polycarbonate plastic). A similar resonance feature 
is found in absorption spectra of O2 (A. Hitchcock, priv. 
comm.). The instrumental features were well determined 
prior to the launch of Chandra in 1999 3 ijPrigozhin et alJ 
1997). 

Along with the known instrumental contribution to the 
oxygen edge, a time-dependent reduction in the effective 
area of ACIS has been detected at low energies. This 
reduction is due to a layer of contaminant that is accu- 
mulating on the ACIS detector which absorbs incoming 
X-rays. The contaminant gradually reduces the detector 
quantum efficiency with time and affects all wavelengths, 
particularly those above 6 A. The Chandra instrument 
teams have studied the spectral profile of the contam- 
inant absorption using LETGS spectra and have mod- 
eled its structure and time-dependence (Marshall ct al. 
2003). An estimated correction for the contaminant is 
now available through the Chandra website 4 . The biggest 
effect is an overall reduction in effective area at the oxy- 
gen edge due to the large optical depth of carbon. Fig- 

3 For in formation on the calibration of the OBF se e 
http: / /www. astro. psu.edu/xray/docs/cal_report/cal_report. html . 
* Kee Tnttp^TT^'ScTrlarvEU'o^eo!!!/^ 



Table 2. Observation Log 



Source Name 


ObsID 


Observation Date 


Time (ks) 


Cyg X-2 


1102 


1999-09-23 


29 


Cyg X-l 


107 


1999-10-19 


9 




3407 


2001-10-28 


21 




3724 


2002-07-30 


14 


4U 1636-53 


105 


1999-10-20 


29 




1939 


2001-03-28 


26 


4U 1735-44 


704 


2000-06-09 


21 


GX 9+9 


703 


2000-08-22 


21 


4U 1543-624 


702 


2000-09-12 


27 


4U 1820-30 


1021 


2001-07-21 


9 




1022 


2001-09-12 


9 



ure 2 shows the change in the effective area over time, 
from 1999 September to 2002 July. In addition, a small 
(r=0.10 in 2002 October) oxygen edge from the con- 
taminant is also detected. There is some discrepancy 
between the predicted optical depth of oxygen (based 
on the behavior of the carbon optical depth) and the 
measured oxygen op tical depth, which is not understood 
l)Marshall et al]l2003T) . Since all of our observations were 
taken before October 2002, the contaminant contribution 
to the oxygen optical depth is <0.10 for all observations, 
which translates to less than 10% of the total instru- 
mental edge. This is comparable to the uncertainty on 
our optical depth measurements (see §5), and therefore 
our results should not be significantly affected by the 
uncertainties in the contaminant contribution. A more 
important issue is the structure of the contaminant edge. 
High-resol ution studies reveal st ructure similar to that of 
the OBF (|Marshall et al.M2003|) . Thus, besides changes 
in the instrumental optical depth, the contaminant also 
adds to the equivalent width of the instrumental absorp- 
tion feature at 23.33 A. It should be noted that the high- 
resolution structure of the oxygen edge was de termined 
using a single observation l|Marshall et al.ll2003|) . and no 
estimate of the error in this measurement is available. 
We assume the contaminant correction is accurate, but 
see §6 for possible problems with the correction. 

In Table 2, we list the observations used in this analy- 
sis. The data sets were reduced using the standard CIAO 
threads 5 . In some cases, the zeroth-order data was not 
telemetered. For these observations, we estimated the 
zeroth-order position by finding the intersection of the 
grating arms and readout streak. After spectral extrac- 
tion, the accuracy of the estimated zeroth-order position 
was verified by comparing the wavelengths of strong in- 
strumental edges in both plus and minus sides of the 
spectra. Events collected during thermonuclear X-ray 
bursts were excluded from this analysis. Since the indi- 
vidual observations of 4U 1820—30 and 4U 1636—53 were 
brief, the spectra from multiple observations were com- 
bined for each source in order to improve statistics. De- 
tector response files (ARFs and RMFs) were created for 
the —1 MEG spectra. We corrected the standard CIAO 
ARFs for the time-dependent change in the response due 
to the contaminant on ACIS using the contamarf tool. 



4. SPECTRAL FITTING PROCEDURE 
5 http:/ /asc. harvard.edu/ciao/threads/ 
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We fit only the 21.5-24.5 A wavelength ranee and 
binned the data to ensure at least 10 counts per bin. An 
absorbed power-law continuum model was used for all 
sources except Cyg X-l, where an absorbed disk black- 
body model was used. We initially allowed the contin- 
uum parameters to vary during the fits, but visual inspec- 
tion revealed that the best-fit parameters underestimated 
the continuum level at wavelengths greater than the edge. 
To better estimate the continuum level, we fit the line- 
free region of the data between 23.7 and 24.7 A and then 
fixed the continuum parameters to the measured values 
in subsequent fits of the oxygen edge region. We will 
present results from the fixed continuum fits, as well as 
note differences between fixed and free continuum meth- 
ods as an estimate of the systematic uncer tainty on our 
result s. The tbvarabs absorption model j Wilms et alJ 
2000) was used with the oxygen abundance set to zero to 
allow for an explicit modeling of the oxygen absorption. 
The equivalent hydrogen column density (Ah) was fixed 
in all fits to a value consiste nt with the measur ed oxy- 
gen edge depth, assuming the I Wilms et al.l 1)2000(1 oxygen 
abundance. 

We modeled the oxygen absorption edge by using two 
edge models which mimic the s tructure of the theo- 
retical determination of the edge l|McLaughlin fc Kirbvl 
119981 iGorczvca fc McLaughlin! 12000(1 . The position of 
the 22.6 A edge was always fixed. For 4U 1735—44 and 
4U 1543—624, the 22.8 A edge position was also fixed in 
the fits. Both edge depths were always varied, although 
for the lower column density sources, only upper limits 
on the depth of the 22.6 A edge were obtained. We al- 
lowed the two edge depths to vary independently in order 
to test if both edge co mponents are detectable as claimed 
bv lTakei et aD {2002ft . Since the ratio of the edge depths 
is a fixed quantity given by the atomic cross-section, we 
should find a consistent value for this ratio if both edges 
are detected in our data. This test allows us to under- 
stand the sensitivity of our data to the edge structure. 

To determine the most appropriate model for the oxy- 
gen edge region, we used the observations of the bright- 
est sources, Cyg X-2 and Cyg X-l (ObsID 107) as tem- 
plates (Figure 3). Our model included absorption lines 
for the 23.5 A and 23. 3 6 A lines seen in other analy- 
ses jPaerels et all l200lt iSchulz et all 120021 iTakei et alJ 
2002). For most sources the widths of these lines were 
allowed to vary along with the line position and depth. In 
some cases, the widths had to be constrained to prevent 
unrealistically large values. The position of the 23.36 A 
line in the fit of the 4U 1543—624 spectra was fixed. 
We also included two absorption lines in the 23.0-23.2 A 
wavelength range and another at 22.8-22.9 A. For these 
lines, we fixed the line width to the instrumental reso- 
lution of 0.023 A FWHM (or a sigma of 0.01 A), while 
allowing the position and depth to find the best-fit val- 
ues. The best-fit parameters were then used to determine 
the equivalent widths (EWs) of the lines. Figures 4 and 5 
show the spectra and best- fit model for all datasets. All 
error estimates are 90%-confidence levels, unless other- 
wise noted, and do not include instrumental errors. The 
absolute wavelength accuracy of the MEG is ±0.011 A 
FWHM. 

5. RESULTS AND ANALYSIS 



By fitting with two edge models, we hoped to re- 
solve the S = 1/2 and S — 3/2 series limits as seen 
in theoretical models. The lack of a consistent ra- 
tio between the two edges, suggests that our data are 
not good enough to consistently detect the S = 1/2 
edge. Since this is true, we have developed a method 
to determine the column density of each edge that 
is independent of the model used to fit the oxygen 
edge. We used the optical depth at 21.7 A, calculated 
from our best-fit model, and the corresponding cross- 
section at this wavelength (Act = 4.2746xl0 -19 cm 2 
where Act is the difference in cross-section from the 
top of the edge to the v alue at A = 21.7 A; 
IGorczvca fe McLau ghlin 200jJ) to determine the column 
density of oxygen (No). We estimate that this method 
has a systematic error of less than 5%. Away from the 
ne ar-edge structure, the ene r gy de pendent cross-section 
of IGorczvca fc McLaughlin! (I2000f) matches the cross- 
section calcu l ated by iVerner fc Yakovlevl (^95). The 
iHenke et al.l (^993) cross-section is larger by rj 10%, 
resulting in small er values of Ao and Ah when the 
IHenke et al.l (^393) cross-section is used for the same op- 
tical depth. It is also important to note that discrep- 
ancies in the inferred column densities may also arise 
when the cross-section at the ionization wavelength is 
used with an edge model placed at an ar bitrary wave- 
length . Using the oxygen ISM abundance of Wilms et alJ 
(2000), we calculated Ah- Table 3 gives the best fit val- 
ues for the edge positions, depths and inferred column 
densities in the fixed continuum models. In addition, 
we give the difference in the calculated Ao between the 
fixed and free continuum models. These column densi- 
ties do not include any correction for oxygen bound in 
dust grains or other molecules. 

Table 4 lists the best-fit values of the line positions, 
widths, and EWs of the oxygen features in the fixed con- 
tinuum method, in addition we show the difference in 
the EW as determined by the fixed and free continuum 
methods. The difference in EW is comparable to the sta- 
tistical error on the measurements. Therefore, our fitting 
method should not greatly affect the results. The 23.51 A 
line is accepted to be the ls-2p transition of neutral oxy- 
gen. By combining the large number of sources in our 
sample, we find that the weighted mean value of the line 
position is 23.508±0.003 A. Fi gure 6 shows the position 
determinations of the O I line for the various datasets. 
Our sample includes one outlier, found in the Cyg X-l 
ObsID 107 spectrum. An examination of the spectrum 
of this source (see Figure 3) suggests that the actual line 
position is higher, consistent with the other sources. The 
mean value of the O I line deviates by only 0.0005 A when 
this observation is excluded from the sample. 

Using our mean value for the O I absorption line as a 
benchmark, we can now shift all of the theoretical data 
to match, allowing us to compare with the other features 
in the oxygen edge region. The weighted mean position 
of the first edge is 22.836±0.007 A. This position is con- 
sistent with the S = 3/2 resonance series predicted at 
22.76-22.84 A. In addition, the lines with best-fit posi- 
tions of 22.86-22.92 A are reasonably consistent with the 
ls-3p resonance line of the oxygen edge at 22.884 A. With 
the resolution of the HETG, we should resolve the ls-3p, 
which must be in the spectra given the detection of the 
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Table 3. Oxygen Edge Parameters 



Source 


A (A) 


T 


N (10 18 cm" 2 ) 


AN (10 18 cm- 2 ) 


Nn (10 21 cm" 2 ) 


Cyg X-2 


22.821±0.013 
22.631 


0.55±0.11 
<0.04 


1.1±0.2 




2.3+0.5 


4U 1636-53 


22.821±0.017 
22.631 


— U. 1 ( 

<0.4 


— U.d 


-0.5 


5.3 _Q <J 




22.631 


U.UJiU.U / 

<0.07 


1 " 31 -0.14 


n ia 


^•'-0.3 


All 1 73^ AA 


99 S^9 

22.631 


i i _i_n a 
0.6+0.4 




— 1 .D 


7-1-9 


GX 9+9 


22.86±0.02 
22.631 


1.1+0.3 
<0.5 


9 9+1. a 

z -0.6 


—0.7 


4 5+ 2 - 4 


4U 1543-624 


22.832 
22.631 


1.4±0.6 
<1.0 




+0.2 




Cyg X-l 
ObsID 107 


22.854±0.017 
22.631 


2.2±0.5 
0.9±0.5 


6.3±1.4 


-0.7 


13+3 


Cyg X-l 
ObsID 3407 


22.838±0.013 
22.631 


0.70±0.10 
0.48+0.11 


2.4±0.3 


+0.1 


4.9+0.6 


Cyg X-l 
ObsID 3724 


99 »qH+ U - U l 7 
ZZ.OOO_q 

22.631 


0.75±0.10 
0.66+0.11 


2.9±0.3 


-0.3 


5.8+0.6 




22.5 23.0 23.5 24.0 

Wavelength (A) 



Fig. 3. — Flux spectra of the oxygen edge region for Cyg X-2 
and Cyg X-l ObsID 107. Line identifications are also shown. 



ls-2p line. We therefore identify the line at 22.88 A as 
the ls-3p resonance line from neutral oxygen absorption 
(see Figure 3 for identifications). 

The combination of line EWs and column densities 
from the edges allows us to study the curve of growth 
of the neutral oxygen lines. A curve of growth describes 
the relationship between the line EW and the absorbing 
column density and can be separated into 3 regimes of 
column density dependence: linear, flat, and square-root. 
In the linear regime, the absorption line is not saturated 



and the EW does not depend on the Doppler broadening 
of the line, yielding a linear dependence on the column 
density (No)- As the line saturates with increasing No, 
the EW becomes explicitly dependent on the Doppler 
broadening as well as on the column density. This is 
known as the "flat part" of the curve of growth, where 
the EW is proportional to (lnAo) 1 / 2 . Finally, as the 
column density increases further, the Lorentzian wings 
of the line profile dominate. In this regime, the EW 
is proportional to (No) 1 / 2 - The flat part of the curve 
of growth is a powerful tool for studying the Doppler 
broadening of the absorbing m aterial- 
Usin g the cross-sections of iGorczvca fc McLaughlin! 
(2000), we have calculated the theoretical curve of 
growth for the ls-2p and ls-3p transitions of neutral oxy- 
gen with velocity dispersions of 20, 100, and 200 km s _1 
(Figure 7). Overlaid are the data from this analysis. 
Both the ls-2p and ls-3p measured EWs are consis- 
tent with the predicted values from the theoretical cross- 
sections, which corroborates our identifications for these 
lines. Our sample of X-ray binaries probes Ao=10 18 - 
10 19 cm~ 2 . For the ls-2p transition (upper panel), this 
places us at the beginning of the (No) 1 / 2 region of the 
curve of growth, where the contribution from Doppler 
broadening is less significant. There seems to be a sys- 
tematic offset of the data from the expected values, which 
is possibly attributable to the use of a Gaussian function 
to fit the lines, instead of a Lorentzian or Voigt profile. 
The ls-3p transition data points are a better match to 
the predicted values (lower panel). For this transition, 
we have moved fully into the flat region of the curve of 
growth. While the error in the data is comparable to the 
given spread in the curve of growth due to differences in 
the Doppler broadening, we can constrain the velocity 
dispersion of the ISM to <200 km s -1 . 

For the line at «23.36 A, our data give a weighted mean 
position of 23.351±0.003 A. Previous authors have at- 
tributed this line to the ls-2p transition of neutral oxygen 
bound in compounds l l Paerels et alJ 120011 iSchulz et alJ 
2002; Tak ei et al J 12002). However, recent calculations 
predict the O II ls-2p absorption line at 23.33 A 
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(jPradhan et al.ll2003j) . and lKawatsura et all (2002) mea- 
sure the O II line at 23.35 A. We therefore attribute this 
line to O II in the ISM (Figure 3). 

Finally, we turn to the two absorption features in the 
23.00-23.15 A range. The weighted mean wavelengths of 
these features are 23.112 A and 23.034 A (with an error 
of 0.004 A). As shown in Table 1, the positions of the 
O III lines are 22.99, 23.08, and 23.14 A. The resolution 
of the MEG is high enough that these lines should be 
resolvable, but with the low countrates, we have binned 
the data 1.5-3 times the instrument resolution, which 
will cause the three lines to blend. Given that we are fit- 
ting three lines with two Gaussian absorption features, 
we would expect the measured positions to lie in between 
the predicted positions, which is what we found. As a 
test case, we fit one of the Cyg X-l spectra (ObsID 3724) 
with three Gaussians for the O III region to determine 
if the lines could be reliably separated. When the wave- 
lengths of the lines were fixed, significant detections of 
absorption lines at those positions were found. When 
the positions were free to vary, the lines were also sig- 
nificantly detected and with wave l engths consistent with 
the predictions of IPradhan et all l|2003f) . although with 
larger errors than found when only two lines were used to 
model the O III absorption. We therefore identify these 
absorption features with O III absorption (Figure 3). 

We used the re sonance oscillator strengths of 
IPradhan et alJ l)2003f) to compute the curve of growth 
for the ionized oxygen lines. For O III, the total EW 
of the lines was used for comparison to the data. A 
curve of growth relates the EW and column density of 
a particular element and ionization state. We have not 
measured the column densities for ionized oxygen, in- 
stead, we compare the EW of the O II and O III lines to 
the column density of neutral oxygen. If we define the 
abundance A = Nx/No, where X is the ionization state 
under consideration, we can relate the column density 
of ionized oxygen to that of neutral oxygen. In Figure 
8, we have plotted the curve of growth for O II (upper 
panel) and O III (lower panel) for three velocity disper- 
sions (20, 100, and 200 km s _1 ) and three abundances 
(A=1.0, 0.1, and 0.001). Both O II and O III are in the 
flat part of the curve of growth, making a determination 
of the abundance more difficult. For reasonable velocity 
dispersions (<200 km s _1 ), we can constrain the abun- 
dance of both ionization states. We find O II/O IwO.l 
andOHI/Ol<0.1. 

It is interesting to note that the data points for O II 
show a large scatter. The two observations with the 
largest deviation are the most recent Cyg X-l observa- 
tions (ObsIDs 3407 and 3724), which are the most af- 
fected by the time-dependent instrumental contamina- 
tion. Therefore, we suggest that the correction for the 
instrumental contamination may not be completely accu- 
rate in its fine structure, and care should be taken when 
drawing conclusions from the most recent data. 

6. DISCUSSION 

We have presented the first global study of oxygen K- 
shell absorption in the ISM. Using HETGS spectra of 
seven X-ray binaries, we have developed an empirical 
model consisting of 2 absorption edges and 5 absorption 
lines. The positions and strengths of the features were 
then compared to theoretical predictions for neutral and 



ionized oxygen features. The neutral oxygen edge was 
resolved into the ls-2p and ls-3p transitions, while two 
edge models were used to describe the remaining, unre- 
solved edge structure. From the curve of growth analysis, 
we constrain the velocity dispersion to <200 km s _1 for 
neutral oxygen in the ISM. We also detect several absorp- 
tion features which we identify with the \s-2p transitions 
of O II and O III. The ISM abundances for O II and O III, 
relative to O I, are found to be «0.1 and <0.1, respec- 
tively. This is the first measurement of the O II/O I and 
O III/O I abundance of the ISM. 

Our analysis relies on using the measured position of 
the neutral oxygen ls-2p line as a benchmark with which 
to shift the theoretical wavelength predictions. Previ- 
ous comparisons of theoretical predictions with labora- 
tory data have shown good agreement of the line po- 
sitions and strengths once a small overall wavelength 
offset (<50 mA) is applied ^McLaughlin fc Kirbvlll99St 
iGorczvca fc McLaughlirJl2000^ . This is at least partially 
due to systematic uncertainties in the absolute wave- 
length calibration of the laboratory data. Laboratory 
measurements of the ls-2p transi tion of neutral ox ygen 
range from 23.489 to 23.536 A fsee lStolte et al.ll997l and 
references therein) with quoted errors much less than the 
discrepancy among the measurements. In contrast, our 
results give a repeatable position for the ls-2p transition 
(over all our sources) due to the HETGS 's well-calibrated 
absolute wavelength scale. With a wavelength accuracy 
of 11 mA, our Chandra data is the best measurement 
of the absolute positions of these features. As an addi- 
tional check, we can compare our resul ts with the labo- 
ratory data for singly ionized oxygen. iKawatsura et all 
(2002) measure the ls-2p transition of O II at 23.35 A, 
0.08 A greater than predicted bv IPradhan et alJ i|2003f) 
and consistent with our measured value. The laboratory 
data supports the need for shifts to the theoretical pre- 
dictions. 

6.1. Previous Results 

iTakei et all l|2002fl fit the spectrum of Cyg X-2 with 
three edges in the oxygen X-shell region, two to represent 
the 2 P and 4 P edges found in atomic oxygen absorption 
and the other for oxygen in molecular compounds. The 
molecular edge had an edge wavelength of 23.13 A, sim- 
ilar to the position found bv lSchulz et alJ l)2002f) for the 
oxygen edge in Cyg X-l. This higher wavelength edge 
coincides with our identified O III absorption features, 
which could resemble an edge featu re in lower resolutio n 
spectra. In the Cyg X-l spectrum of lSchulz eFaD(|2Q02), 
there is a significant count excess on the shorter wave- 
length side of the edge and at the time, this excess was 
attributed to either an instrumental effect or to emission 
lines from the source. Our model is more successful in 
fitting this struct ure (see Figure 3). 

Other authors jPaerels et al.ll200lt iSchulz et al1l2002t 
ITakei et al.ll2002|) have previously identified the O I ls- 
2p line and the line at 23.35 A in the spectra of X-ray 
binaries. The 23.35 A line has been attributed to the 
ls-2p transition in iron oxides based on the results of 
IWu et all l)1997j) which suggest a narrow absorption fea- 
ture at «23.3 A. But in addition to the narrow absorp- 
tion feature, th ere is a broad edge feature at «22.9 A 
tWn et alJH997l) . In fact, given the relative cross-section 
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Fig. 4. — Flux spectra for six X-ray binaries included in this study. The best- fit model is also shown. The data have been smoothed for 
illustrative purposes and representative error bars are shown. The dashed lines indicate the positions of the identified features. 



between the narrow line and the edge, the edge should be 
detectable before the line. This is opposite of the atomic 
oxygen cross-sections, where the narrow ls-2p lines have 
larger cross-sections than the associated edges. This edge 
feature is not required by our data, which is fit well by 
just the neutral, atomic oxygen edge. Finally, the struc- 
ture of molecular edges is dependent on the com position 
of the dust in the ISM, which is highly uncertain {brainc 
2003). This, coupled with a lack of laboratory and theo- 
retical results for molecular oxygen cross-sections, makes 
predicting the position and strength of molecular line 
features difficult. 

We therefore attribute the 23.35 A feature to ab- 
sorption from O II since we expect ionized oxygen to 
be present in the ISM at the position we determined. 
In the warm, ionized phase of the ISM, the ioniza- 
tion of oxygen should track th at of hydrogen due to a 
strong charge exchange process (iField fc Steigmanll97lD . 



Therefore we wo uld expect t h at N (Q II)/iV(0 I) « 
N(H Il)/iV(H I). iTakei et alJ (12002ft estimate th at for 
CygX-2, N(R Il)/iV(H I) = 0.32. IRevnoldl iTlQQI found 
H II/H I of 0.2-0.6 along different lines-of-sight through 
the galaxy. Our estimate for N(0 Il)/iV(0 I) w 0.1 is less 
than that found for hydrogen, but due to multiple elec- 
tron ejections during photoionization (|Watson fc Kunzl 
1975?) , a more accurate comparison should be [N(0 II) + 
N(0 III)]/AT(0 I). In either case, the presence of O II is 
expected at roughly the measured value, supporting our 
identification. 

6.2. Contributions from Molecules and Dust? 

Our data are well fit by our model which includes only 
neutral and ionized oxygen contributions. But the ISM 
should also contain oxygen in various compounds, includ- 
ing carbon monoxide (CO) and dust. In X-ray spectra, 
molecular forms of oxygen should produce a \s-2p reso- 
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Table 4. Oxygen Line Parameters 



Source 


A (A) 


o- (A) 


EW (A) 


AEW (A) 


Off X-2 


23 508+0 004 


026+0 005 


067+0 01 1 






Zo.OOZlCU.UUO 


U.U10IEU.UU t 


n nAi -t-n m 1 
u.u^izcu.ui 1 






9Q r\"7Q-i-C\ ni o 

zo.u / oztu.uiy 


n ni 
U.U1 


n nno4-n nns 
U.UUyztU.UUo 






zo.uoinu.uo 


U.Ul 


^n m 7 






22.899t .oi7 


0.01 


0.019+0.011 




4U 1636-53 


23.507+0.011 
23.340+0.008 


040 +u ' ui;1 
u.utu_ 014 

<0.02 


0.09+0.03 
0.031+0.017 


-0.03 
-0.003 




23.102+00!. 4 

—0.06 


0.01 


0.020+0.010 


-0.005 




23.056+0.013 


0.01 


0.027+0.010 


-0.004 




22.915+0.013 


0.01 


0.034+0.010 


-0.004 


4U 1820-30 


23.514+0.010 


0.034+0.015 


0.07+0.02 


-0.02 




23.351+0.014 


<0.04 


0.04+0.02 


-0.01 




23.086+0.017 


0.01 


0.033+0.013 


-0.004 




23.03+0.06 


0.01 


0.023+0.013 


-0.004 




22.867+0.015 


0.01 


0.036+0.016 


-0.004 


4U 1735-44 


23.503+0.009 


0.033+0.015 


0.08+0.03 


-0.02 




23.354+0.017 


0.027+0.014 


0.06+0.02 


-0.03 




23.095+0.017 


0.01 


0.027+0.013 


-0.016 




23.020+0.007 


0.01 


0.036+0.009 


-0.010 




22.861+0.006 


0.01 


0.040+0.009 


-0.008 


GX 9+9 


23.517+0.009 


0.036+0.016 


0.08+0.03 


-0.02 




23.367+0.017 


0.026+0.014 


0.050+0.019 


-0.018 




23.133+0.010 


0.01 


0.024+0.010 


-0.007 




23.017+0.011 


0.01 


0.034+0.010 


-0.006 




22.906+0.018 


0.01 


0.031+0.011 


-0.006 


4U 1543—624 


23.507+0.011 


0.025 


0.068+0.011 


-0.005 




23.33 


0.01 


0.026+0.017 


-0.009 




Zo.14iliU.UlU 


U.Ul 


n nQe;4-n m ^ 
U.UoDztU.Ulo 


n m 
-U.Ulz 




oq n^n^n ni 1 
zo.UoU±U.Ull 


n ni 
U.Ul 


n nQC4-n m ^ 
U.UooztU.Ulo 


ri (11(1 
-U.U1U 




00 C7i fine 
zz.o 1 l±U.UUo 


n ni 
U.Ul 


U.U4o±U.U10 


n nno 
-u.uuy 


Cyg X-l 


23. 482±0. 019 


0.086_ 002 


U.21±0.0o 


-0.06 


ObsID 107 


23.333+0.015 


<0.05 


n5+ 05 


-0.001 




23.106+0.017 


0.01 


0.044+0.011 


-0.009 




23.036+0.009 


0.01 


0.056+0.006 


-0.002 




22.910+0.013 


0.01 


0.053+0.011 


-0.008 




OQ £1 QJ-d nno 


n n^i -4-n ni 7 


n nQ^-i-n m 7 


-U.Ulo 


ObsID 3407 


23.370+0.012 


019+ 018 
u.ui»_ on 


0.020+0.011 


-0.001 




23.092+0.009 


0.01 


0.018+0.005 


-0.001 




23.044+0.010 


0.01 


0.017+0.006 


-0.002 




22.885+0.017 


0.01 


0.011+0.006 


-0.002 


Cyg X-l 


23.505+0.006 


0.05 


0.094+0.007 


-0.005 


ObsID 3724 


23.349+0.010 


0.01 


0.009+0.005 


-0.007 




23.118+0.010 


0.01 


0.030+0.006 


-0.007 




23.039+0.018 


0.01 


0.024+0.006 


-0.007 




22.870t°;g} 8 


0.01 


0.022+0.014 


-0.006 



nance feature at a wavelength shifted from that of atomic 
oxygen due to the effects of the bonds on the electron en- 
ergy levels. We recognize that there are unidentified fea- 
tures in the spectra with strengths comparable to some 
of the smallest EW features we identify. It is possible 
that these features could be attributable to molecular 
features, though any quantitative statements depend on 
predictions or measurements of the strength and posi- 
tions of these features, which are unavailable at this time. 

Although our data do not require such contributions, 
we can estimate an upper limit to the amount allowed by 
our data. If we take the error estimate for our measured 
column densities to be the upper limit on the amount 
of oxygen bound in molecules, and assuming the cross- 
section for molecular oxygen is comparable to that for 
atomic, then 10-40% of the oxygen can be bound and 
still not be detected. A similar limit can be placed if we 
identify the 23.35 A feature as due to a molecular con- 
tribution as has been previously suggested ijPaerels et alJ 



l200lt ISchulz et al.ll200l iTakei et alJl200^ . In this case, 
the ratio of molecular to atomic oxygen would be roughly 
the same as we suggest for O II/O I, «10%. This as- 
sumes that the oscillator strength for molecular oxygen, 
which is unknown, is similar to that of O II. These es- 
timates compa re well with t he ex pected amount from 
other studies. lAndre et alJ l)2003j) calculated that the 
theoretical maximum amount of oxygen which can be 
found in dust grains is (O/H)^^ < 180 ppm by com- 
paring the ISM abundance of oxygen to magnesium, sil- 
icon, and iron. Given the measured oxygen abundances 
in the ISM and the Sun, they found an average for 
(0/H)dust = 109 p pm, which translat es to 27% of the 
oxygen m the ISM. lAndre et all l)2003j) also found that 
the amount of oxygen bound in grains varied between 
lines of sight, from 1-40%. Higher signal-to-noise spec- 
tra will be required to measure the effect of dust and 
molecules on the oxygen edge in X-ray spectra. 
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Fig. 5. — Flux spectra and best-fit model for the three Cyg 
X-l observations. The data have been smoothed for illustrative 
purposes and representative error bars are shown. The dashed lines 
indicate the positions of the identified features. Note the reduction 
in the relative strength of the O II to O I absorption lines (23.35 
and 23. 5i A, respectively). We attribute this to possible problems 
in the modeling of the time-dependent contaminant. 



6.3. Alternatives to Interstellar Origin 

Although we have attributed the identified features to 
interstellar oxygen absorption, we also consider alterna- 
tive explanations for completeness. Pre-launch studies 
of the instrumental oxygen absorption features do not 
show any significant features that would account for the 
absorption lines we have found, other than the 23.33 A 
absorption line from the OBF ijPrigozhin et al.l fl997'). 
In addition, the correlation of the line EWs with neutral 
column density would not be seen from an instrumen- 
tal feature. Instrumental features would be expected to 
give the same EW regardless of the source and would be 
seen in observations of all sources, including low column 
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Fig. 6. — Best fit position of the ls-2p absorption line from O I. 
The weighted mean value of the line position is given. The outlier is 
from the Cyg X-i ObsID 107 data. All quoted errors are statistical 
errors only and do not include the instrumental wavelength error 
of ±0.011 A (FWHM). 
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Fig. 7. — Curve of growth for the ls-2p and ls-3p transitions of 
neutral oxygen. The solid, dashed, and dotted lines indicate the 
predicted EWs for velocity dispersions of 20, 100, and 200 km s —1 , 
respectively. 
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density extragalactic sources, like those used for calibra- 
tion of the contaminant contribution. We know of no 
identification of these features in lower column sources. 
Therefore, we reject an instrumental origin for our de- 
tected line features. 

That being said, there does seem to be some effect 
on the EW determinations of the O II line from the 
time-dependent contaminant absorption. For most of 
our data, the contaminant effect is small and does not 
impact our results for the ISM. The most recent obser- 
vations however, which would have the largest contam- 
inant correction, have lower EWs than predicted from 
the curve of growth of the other data points. One pos- 
sible explanation is that there is a non-negligible O II 
absorption line in the spectra of the calibration source. 
We note that the structure of the oxygen contaminant 
absorption is primari ly taken from a single observation 
({Marshall et al]l2003) . By not accounting for ISM O II, 
the contaminant correction is larger than necessary, re- 
ducing our measured EWs, particularly as the correction 
grows larger. It is hoped that our study can help the 
calibration team understand ISM contributions to the 
oxygen edge region in order to more accurately model 
the contaminant. Additionally, the most recent obser- 
vations of Cyg X-l also show an excess of counts in the 
22.2-22.6 A range. Since these are the only observations 
with this feature, it is again possible that it is due to the 
effect of the contaminant or the contaminant model. 

Accepting that the line features are from ionized oxy- 
gen species, we must ask the question are these features 
truly from the ISM, or rather do they originate in ma- 
terial local to the X-ray binaries. The continuum emis- 
sion from the binaries provide a great source of ionizing 
radiation, which could create ionization bubbles around 
the systems. The amount of ionized oxygen found in 
each system would depend on the luminosity of the X- 
ray source, and the density and composition of material 
surrounding the source. To first order, if the density of 
material is taken to be similar for all systems, the more 
luminous systems will be able to ionize more oxygen and 
would then be expected to have the largest measured 
EWs. Of our sources, Cyg X-l and Cyg X-2 are the 
most luminous, yet Cyg X-2 has the some of the smallest 
measured EWs. This is expected if the ionized mate- 
rial is interstellar, since Cyg X-2 also has the smallest 
neutral oxygen column density from the edge determina- 
tion. While the O II data points do not show a significant 
trend with column density, most likely due problems with 
modeling of the contaminant as mentioned previously, 
the EWs of O III are correlated with the neutral oxy- 
gen column density. Given that the sources included in 
our sample span more than two orders of magnitude in 
luminosity and several orders of magnitude in orbital sep- 
aration, it is highly unlikely that the correlation between 
the O III EW and the neutral column density could be 
produced by absorption from material local to the bina- 
ries. On the other hand, the correlation is exactly what is 
expected if all of the absorption features are interstellar 
in origin. 

Along with the lower ionization species found here, 
we would expect the binaries to produce absorption fea- 
tures from more highly ionized forms of oxygen. Of our 
source s, only Cyg X-l has shown other absorption fea- 
tures llSchulz et al.ll2002t lF7u7et a,U 12003 IMiller et alJ 
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Fig. 8. — Curve of growth for the O II and O III absorption fea- 
tures. The O III EW is the total EW from the multiple absorption 
lines. The solid, dashed, and dotted lines indicate the predicted 
EWs for velocity dispersions of 20, 100, and 200 km s , respec- 
tively. Also included are predicted EWs for various abundances 
(A = N-^/Nq) relative to neutral oxygen. 



2003), and it is the only source with multiple observa- 
tions that can be used to study variations of the column 
densities. In the three observations of Cyg X-l, we mea- 
sure a change of the neutral oxygen column density of 
a factor of 2, which seems unreasonably high for varia- 
tions in the ISM. Cyg X-l is a high-mass X-ray binary 
system that accretes from the stellar wind of the young, 
massive companion. This, combined with the detection 
of other more highly ionized absorption lines, points to 
local material in this system. If the fluxes of the obser- 
vations are compared, we find that ObsID 107 has the 
lowest continuum flux at the oxygen edge region and the 
highest neutral column density. The fluxes are larger in 
ObsIDs 3407 and 3724 with smaller columns. We pos- 
tulate that as the flux of Cyg X-l incre ases, so does the 
ionization of the local material (but see lWen et aLlll99ft 
iFeng et all 12003 for a discussion of the dependence on 
orbital phase). In ObsIDs 3407 and 3724, the local ma- 
terial is most likely highly ionized, removing any local 
contributions to the neutral and low ionization features. 
Measurements of the other absorption lines may help to 
verify this hypothesis. 

This paper represents the first in a series of studies 
of the properties of the ISM using high-resolution X-ray 
spectra. Future work will include studies of lower column 
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density sources to map out the curve of growth of the dif- 
ferent oxygen absorption lines. More data will allow us 
to better constrain the velocity dispersion and ionization 
fraction of the oxygen in the ISM. Additionally, we will 
perform a similar analysis on the other astrophysically 
abundant elements, which can be used to study abun- 
dances and depletion in the ISM. Studies of the ISM in 
X-rays provides a unique opportunity to measure a large 
variety of elements and ionization states. 

Some interesting questions that will be addressed in 
the future are the relationship between the warm, neu- 
tral (WNM) and warm, ionized (WIM) phases of the 
ISM and the source of the ionizing flux. From an ultra- 
violet observation of t he Ga lactic halo star HD 93521, 
iSpitzer fc Fitzpatrickl (1993) found tha t the WIM coin- 
cided with the WNM. Ot her studies <|Howk fc Savage! 
1999; [Reynolds et alJll995|) have found similar results. 
It has been suggested that the ultraviolet flux from O 
stars would provide enough energy to produce the WIM 



(see iHeiles fc Kulkarnl 11987 ). However, the ionizat ion 



fraction of helium found by Reynolds fc Tuft el l|1995() is 
much lower than expected for the O star ionization mod- 
els. New observations from the Far Ultraviolet Spectro- 
scopic Explorer and the Wisconsin H- Alpha Mapper will 
help to resolve this issue. In addition, the high-resolution 
X-ray spectra available from Chandra and XMM will pro- 
vide information on the ionization fraction of heavier el- 
ements in the ISM, providing more data from which to 
understand the source of the ionizing radiation. 
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